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Abstract 

The exact expression for the phase-dependent linear conductance of a weakly 
damped superconducting quantum point contact is obtained. The calculation 
is performed by summing up the complete perturbative series in the coupling 
between the electrodes. The failure of any finite order perturbative expansion 
in the limit of small voltage and small quasi-particle damping is analyzed in 
detail. In the low transmission regime this nonperturbative calculation yields 
a result which is at variance with standard tunnel theory. Our result predicts 
the correct sign of the quasi-particle pair interference term and exhibits an 
unusual phase-dependence at low temperatures in qualitative agreement with 
the available experimental data. 
PACS numbers: 74.50.+r, 85.25.Cp, 73.20.Dx 



Since the early stages in the theory of Josephson junctions it has been customary to 
write the total current through the junction as 

/ = ij sin + C7 o (l + ecos0)V, (1) 

where the total superconducting phase difference between the electrodes is related to the 
applied bias voltage by dfyfdr = 2eV/h = u . 
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Eq. (1) was first derived from a microscopic model by Josephson An equivalent 
expression has been widely used to describe superconducting point contacts 0. The second 
term in Eq. (1) defines a phase- dependent conductance, G(<f)) = G (l + e cos (ft), whose 
existence was confirmed by a series of experiments during the seventies 0. However, it soon 
became apparent that the experimental results for low temperatures seemed to be best fitted 
by Eq. (1) with a value of the parameter e ~ —1, while tunnel theory predicts e = +1 in the 
limit T, V — > [[| . Furthermore, an experiment measuring the complete phase dependence 
of the linear conductance in a point contact showed a strong departure of G(4>) from the 
simple cos^-like form |J, given by tunnel theory. There was at the time a large number of 
theoretical works trying to explain the discrepancy between tunnel theory and experiments, 
most of them relying on the introduction of a phenomenological broadening of the Riedel 
peak (for a review see ref. Q). However, a completely satisfactory explanation of this issue 
seems still to be lacking. 

On the other hand, there has been in the last few years a renewed interest on the theory 
of superconducting weak links associated with the increasing technological capability for 
the fabrication of nano-scale superconducting devices. This opens the possibility of a closer 
comparison between theoretical predictions and clean experiments even on a nearly atomic 
scale @. It therefore seems appropriate at this point to perform a careful re-examination 
of the transport properties in superconducting point contacts going beyond the limits of 
standard tunnel theory. In this direction, one should mention some important theoretical 
contributions JF|||, clarifying the crucial role played by multiple Andreev-reflections in the 
explanation of the sub-harmonic structure in superconducting point contacts. This sub-gap 
structure can be currently analyzed experimentally with increasing resolution |6],|9|]. 

As will be shown in this letter, the presence of edge singularities in the spectral densities 
of the superconductors leads to an enhanced weight of these multiple scattering processes. 
As a consequence, a nonperturbative calculation is needed to obtain the correct result for 
the phase-dependent conductance. In a non-biased junction, the summation of the infinite 
series of multiple scattering events leads to the appearance of bound states inside the super- 
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conducting gap, whose existence has been discussed in previous theoretical works |T0|-|T2 
In the same way, the presence of bound states will allow us to obtain an exact expression 
for G{4>) valid in the limit of small quasiparticle damping. 

We consider for simplicity the case of a symmetrical contact. Then, with L and R 
representing the left and right electrodes, we have |Ax| = \Ar\ = A and <fi = 4>l — 4>r, where 
A is the modulus of the superconducting order parameter, representing the total phase 
difference which is supposed to drop abruptly at the interface between both electrodes. 

To describe the biased point contact we use the following model Hamiltonian 

H = H L + H R + H LR + ^- (ft L - N R ) , (2) 

where Hl and H R are the BCS Hamiltonians for the uncoupled electrodes, (N L ,N R ) being 
the corresponding total number operators, and eV is the applied bias voltage. The term 
Hlr coupling both electrodes is assumed to have the following form 

HlR = y](tii c irr c i<r + tji c ja c ia)> (3) 

where stand for orbitals on the (left, right) electrodes respectively. For the present 

analysis it is convenient to consider first the simplest case in which there is a single channel 
connecting both electrodes (the multi-channel case will be briefly discussed at the end of the 
letter). In this case (i = L,j = R) denote the two orbitals connected by the single hopping 
element Ilr = t. We perform the standard unitary transformation |T5|Jl6|] , by means of 
which the system dynamics is governed by the following time-dependent Hamiltonian 

H(t) = H l + H r + J2(te i<t>iT)/2 c{ a c R<7 + te-^c^CL,,), (4) 
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where 0(r) = <fio + 2eVr/h. Within this representation the time-dependent phase only 
appears in the phase factors multiplying the hopping elements. The transport properties of 
this system can be analyzed using nonequilibrium Green functions techniques P,0,ITB[ with 



the time-dependent coupling term treated as a perturbation. The most relevant quantity in 
this formalism is the nonequilibrium distribution function G + ~, which in a superconducting 
broken symmetry (Nambu) representation is defined by 



/ 



< cU(T')CiAr) > < Co|(V)c,;t(t) > 



\ 



<cI(t')c]At) > <c,At')cI(t) > 



(5) 



In terms of these functions the current through the contact can be written as 



Ht) = X [i(r)G + RL (r,r)-i^r)Gi R (r,r) 
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(6) 



where t is the matrix hopping element in the Nambu representation 



te^( T )/ 2 



V 



o 





-te~ irt>[T)/2 



(7) 



Within this perturbative approach the tunnel theory expression for the current (Eq. (1)), 
can be obtained at the lowest-order in Hlr- The conductance Gq{1 + ecos(/>) thus obtained 
becomes a divergent quantity in the limit V — > [|]]. In order to ensure the existence of a 
linear regime, a finite energy relaxation rate rj must be introduced into this superconducting 
mean field theory (rj represents the damping of the quasi-particle states, which in a real 
system is always present due to inelastic scattering processes). As we shall see, according to 



the value of r\ and the normal transmission coefficient of the junction, a ||13|| , two different 
regimes can be identified: the weakly damped regime, for which rj aA and the strongly 
damped case, where rj ^> aA. In this work we are mostly concerned with the analysis of 
the first regime, where the most interesting effects appear. 

A remarkable fact about the perturbative expansion in the weakly damped situation 
is that contributions corresponding to higher order processes turn out to be increasingly 
divergent in the zero bias limit |I4 |. In particular, it can be easily demonstrated that 
contributions to the total current of order t 2n , n > 2, diverge like ~ t 2n /r] n ~ 1 (the lowest 
order contribution diverges as ~ t 2 In rj) . This result is a direct consequence of the increasing 
contribution from the superconducting gap edges singularities. Therefore, a correct answer 
cannot be found in principle by means of a finite order perturbative expansion. 

One could draw a formal analogy with the case of a high-density electron gas, where 
the diagrammatic expansion in the bare Coulomb potential is also increasingly divergent. 



As in that case, the solution can be found by "dressing" the perturbative potential, i.e. 
Hlr- In the present problem the dressed quantities (left-right coupling, propagators) can be 
exactly obtained in the zero voltage limit by evaluating the complete perturbative series. To 
this end, we find it convenient to express all quantities in terms of a renormalized left-right 
coupling element which satisfies the following Dyson equation 

f ^(r, r') = t(r)6(r - r') + t(r)g a /(r - n )& (n) g a / (n - r 2 )T^(r 2 , r'), (8) 

where gY and gY represent the (advanced, retarded) Green functions of the uncoupled 
left and right electrodes respectively (integration over internal times is implicitly assumed). 
From Eq. (8), the relation between the renormalized coupling T and the exact (advanced 
and retarded) Green functions is easy to obtain. In the same way, the nonequilibrium 
distribution function G + ~ , which is related to G r and G a , can be written in terms of T [JT9j . 



Integral equations like Eq. (8) adopt a simpler form when Fourier transformed with 
respect to their temporal arguments PJI5fl. Defining the Fourier components T n ^ m {uj) as 

f n , m (uj) = fdT f dT > e -^(T)-rnrtT>))/2 e -iu,(T-T')f^ T ,^ (g) 



the total current can then be expressed in the form J(r) = Y, m Im ex P iw4>if) /2, where the 
complex coefficients, J m , do not depend on 0(r) and are given by 

2e 



I 



m 



{ J &" E + n^)T:}jl{u + m^) - g r Mn,jr^ + n ^)tf 
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(10) 



It can be seen from Eq. (8) that T n ^ m = for even n — m and therefore only even Fourier 
components of the current are different from zero. 

For the following analysis it is useful to divide the total current into dissi- 
pative and nondissipative contributions. The supercurrent part, given by 1$ = 
— 2 E m >o I m (-^m) sin[m0(r)], tends to a finite value in the limit V — > 0. On the other 
hand, the dissipative part is given by Id = Iq + 2 J2m>o R- e (-^m) cos[m0(r)], and goes to zero 



as In ~ G(4>)V, G((f>) being the zero voltage conductance. The linear term can be straight- 
forwardly derived from Eq. (10) by expanding the Fermi functions appearing in g~^~ R JET] up 
to first order in V and evaluating the remaining factors at zero voltage. 

In this limit the Fourier components satisfy T n ,m = 2o jm _ n = T m _ n , and can be shown 
to obey the simple recursive relations 

T n+2 {uj) = z(w)f n (u) 
f_ n _ 2 (o;) = z{u)f_ n (u) (n > 1), (11) 

where z(u) is a scalar complex function. In the weakly damped regime and within the energy 
interval A > \u\ > Ay/1 — a this function reduces to a phase factor z(u) = expi<f(u), where 
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(p{u) = arcsin ^—y/A 2 -u 2 ^u 2 - (1 - a)A 2 j . (12) 

This clearly shows that in the weakly damped regime and within this energy interval all 
multiple scattering processes become equally important. Therefore, all Fourier components 
contribute to the renormalized coupling in this region, giving rise to singularities which 
can be shown to be associated with the existence of interface bound states. In fact, the 
renormalized coupling in this energy region can be easily obtained from Eqs. (11) and (12), 
giving 

VT n ( W )e ! ^ 2 = Il[LJ)€ + * , (13) 

which exhibits singularities at ip(uj) = ±0. From Eq. (12) it follows that these singularities 



correspond to simple poles at us = ±Ay 1 — asin 2 (0/2). These are the interface bound 
states inside the gap of a superconducting point contact, as derived by different authors 

EMU- 



In the same way, the complete harmonic series must be evaluated in order to obtain the 
contributions to both the dissipative and nondissipative parts of the current coming from the 
energy range A > |o;| > Ay/1 — a. Again, these infinite summations can be easily performed 
making use of the recursive relations of Eq. (11). It is then found that the integrand for 



both parts of the current becomes singular at lu = ±ujs- The contribution of these poles 
yields 

h{4>) = ijr i = =tanh(^) (14) 

2h yfl- a sin 2 (0/2) 2 

and 



1d{4>) 



h 16/7 



Act sin ,/3ujs* 
=sech(— — 



l-asin 2 (0/2) 
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(3V. (15) 



In Eq. (14) the previously known result for the zero bias supercurrent is recovered |10|-[T2 



The expression for the dissipative current given above is the main result of this letter. The 
linear conductance thus obtained can be seen to depend on r\ as ~ I/77, i.e. proportional to 
a relaxation time. This result seems more physically sound than that of tunnel theory which 
predicts a dependence ~ In 77. Notice that in the low barrier transparency regime Eq. (15) 
depends on the transmission coefficient as a 2 which means that Andreev reflection processes 
dominates over single quasiparticle tunneling in the zero voltage limit. 

It is also worth commenting that Eq. (15) can be derived in a different way, by relating 
the linear conductance to the equilibrium current fluctuations via the fluctuation- dissipation 
theorem, giving further support for the validity of this expression. This will be the subject 
of a forthcoming publication. 

Our theory yields a phase-dependent linear conductance which strongly deviates from the 
tunnel theory result of Eq. (1). In the limit of low barrier transparency, Eq. (15) predicts 
G(4>) ~ 1 — cos(20) instead of G(<f>) ~ 1 + ecos</> of standard tunnel theory. Therefore, the 
linear conductance in Eq. (1) can never be recovered in the weakly damped regime. On 
the other hand, with increasing values of 77 multiple scattering processes are progressively 
damped (the function z(u) is no longer a phase factor decaying exponentially with 77); 
eventually, when rj 3> aA only the lowest-order processes contribute to the current and Eq. 
(1) is recovered. This explains the discrepancy between tunnel theory and the experiments 
(usually referred to as the "cos0 problem"), because the experimental conditions should 
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correspond to the weakly damped case in order that the Josephson effects could be observed 

Another interesting limiting case of Eq. (15) corresponds to the ballistic, i.e. a — > 1, 
regime. In this case and for large temperatures G(<p) behaves approximately as (1 — cos0), 
in agreement with the result given by Zaitsev [Q. However, the most unusual phase- 
dependence of G((p) appears for high values of the transmission and low temperatures (ksT < 
A). This is illustrated in Fig. 1, where G(<f)) is plotted for two different temperatures and 
increasing values of the transmission. The only experiment where the full phase dependence 
of G(4>) was measured is, to our knowledge, that of ref. [f|. Their measured G(<f)) strongly 
deviates from a cos 0-like form, being almost negligible for small values of <j) and exhibiting a 
large increase around ~ 7r/2. As can be observed in Fig. 1, this behavior is in qualitative 
agreement with our results at any given temperature for sufficiently large transmission. 
However, a detailed comparison should require a more exhaustive experimental study of 
G(4>) for different barrier transparencies and temperature regimes. We believe that these 
measurements are now becoming feasible with recent advances in the fabrication of nanoscale 
superconducting contacts. 

The multi-channel generalization of our results is formally straightforward. For a general 
contact geometry it would lead to a superposition of contributions like those of Eqs. (14) 
and (15) for each transverse mode, which can in principle have different transmission prob- 
abilities. We do not expect that this effective averaging process would alter in a significant 
way the phase dependence of Eqs. (14) and (15). 

In conclusion, it has been shown that a nonperturbative calculation is needed for ob- 
taining the total current through a weakly damped superconducting point contact in the 
linear regime. Using a simple model Hamiltonian we are able to obtain exactly the phase- 
dependent linear conductance. The resulting expression is in good agreement with the 
available experimental data and we believe it can provide a motivation for more detailed 
experimental studies. 
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FIGURES 

FIG. 1. Phase dependence of the linear conductance given by Eq. (15) for two different tem- 
peratures and increasing values of the normal transmission coefficient (G((p) is normalized to its 
maximum value). 
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